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ABSTRACT 

J2 f 3 ?? 

Some theoretical and experimental developments of the last decade are 
discussed with emphasis on techniques that can be used to calculate viscosity, 
thermal conductivity, and diffusion coefficients assuming appropriate interatomic 
of intermolecular force laws. Three main topics are considered: First, inter- 

molecular potential energy functions for which collision integrals are now avail- 
able; second, heat conduction in chemically reacting gases; and finally, a closely 
related topic, the heat conductivity of polyatomic gases. ftcA&fosx-- 

INTRODUCTION 

This paper considers theoretical and experimental developments concerning 
the transport properties since about 1954; in other words, since the publication 
of "Molecular Theory of Gases and Liquids" by Hirschf elder, Curtiss, and Bird' 1 ', 
which effectively summarizes most prior work. Furthermore, the discussion is 
limited to techniques that have engineering usefulness - techniques that can be 
used to actually compute transport properties algebraically. 

Thus, three main topics are discussed: First, the collision integrals, 

basic to transport property calculations, that are now available for a wide 
variety of intermolecular force laws; secondly, heat conduction in chemically 
reacting gases followed; finally, by consideration of a closely related subject, 
the thermal conductivity of polyatomic gases. 

In general this review will make reference to the literature rather than 
presenting computational methods as such. Equations will be avoided except 


insofar as they serve to illustrate tl 
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in order to carry out any actual calculations, the reader must consult the 
literature. 

POTENTIAL ENERGY FUNCTIONS AND THE COLLISION INTEGRALS 
The rigorous Chapman-Enskog theory yields the following expressions for the 
transport properties of dilute monatomic gases^: 

Viscosity 


Thermal conductivity 


Self -diffusion coefficient 
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These formulas involve quantities such as the atomic mass m, the Boltzmann 

constant k, the temperature T, the heat capacity c y (= ^ k), and the density 

p, which are well known. However, the formulas also contain cross sections, or 
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more properly collision integrals a fl' * ' and 0 (! l ' , and to compute 

these the intermolecular force law must be known. 

For a spherically symmetric potential, which may be written in a dimension- 
less form as 


ikl = cp* = f j = f (r*) (4) 

the collision integrals are obtained by a triple integration. (Here e is an 
energy and a is a distance characteristic of the potential. ) First it is 
necessary to compute the angle of deflection: 
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x(g*,b*) = 
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where b* = b/a is the reduced impact parameter. (The impact parameter b is 
the distance of closest approach in the absence of the potential cp. ) Further, 
r* = r / a, where r m is the distance of closest approach in the presence of 
the potential, and g* 2 = i mg 2 / e is the reduced relative kinetic energy (g 
is the initial relative speed of the colliding molecules). 

Once the angle of deflection has been obtained as a function of g* and b*, 
a velocity-dependent cross section is computed: 


Q (z) *(g)* = 
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(1 - cos z X)b*db* 
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Finally, the Q 
weighting factor: 
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are averaged over all velocities, with an appropriate 
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Thus the collision integrals are a function of reduced temperature 


T* 5 kT/e. 


The purpose of presenting Eqs. (5) - (7) is to show the manner in which the 
potential energy cp influences the transport properties. The potential appears 
explicitly in the integrand for the angle of deflection and is then averaged by 
three integrations. Consequently, the collision integrals are insensitive to 
the details of the intermolecular potential, and we cannot expect experimental 
transport property measurements to provide much information about the force law. 
To put it another way, the experimental data cannot be used to determine the 
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potential; the data can only he used to veto candidate potentials. 

For most potential energy functions Eqs. (5) - (7) must he evaluated numeri- 
cally. The first such calculations for a realistic force law - the Lennard- Jones 
(12-6) potential - were carried out independently hy four different groups 0 . 
Furthermore, the Wisconsin groups showed that experimental viscosity data could 
he fit rather well to yield molecular constants o and € in reasonable agree- 
ment with corresponding parameters obtained from equation of state data. 

The Lennard-Jones (12-6) potential 

<P(r) - 4€ [(?f - (?f] (8) 

combines an inverse sixth power attractive potential with an inverse twelfth 
power repulsion. The attractive portion has a theoretical basis in the dispersion 
forces, hut the twelfth power repulsion was chosen merely for mathematical con- 
venience. In view of the considerable success of the Lennard-Jones (12-6) 
potential it is not surprising that the next step was to build some flexibility 
into the repulsive part of the potential. This was accomplished by introducing 
an exponential repulsion in place of the inverse twelfth power repulsion; the 
collision integrals for these potentials, shown in Fig. 1, have been computed 
by Mason * . The Lennard-Jones (12-6) potential is also shown in Fig. 1 for 
purposes of comparison. 

The exponential-6 potential serves very well for the noble gases and other 
molecules that are approximately spherical. Figure 2 shows experimental viscosity 
data for argon over a wide temperature range; a curve computed for rigid elastic 
spheres is shown as a dashed line. The potential functions are also shown in 

D 

Fig. 2. The exponential-6 potential has been chosen to take account of equa- 
tion of state and crystal properties as well as viscosity coefficients. For 
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example, the interatamle distance in solid argon is 3.8 Angstroms - just inside 
of the minimum of the potential energy curve. The same potential also provides 
a satisfactory explanation of the thermal conductivity (Fig. 3) and self- 
diffusion coefficient (Fig. 4) for argon. 

Since theoretical calculations have been very successful in describing the 
transport properties for valence- saturated nonpolar gases, it is not surprising 
that the collision integrals have now been computed for potentials applicable to 
other situations. 

Molecules possessing permanent dipole moments have an intermolecular potential 
that is not spherically symmetric. The dipoles give rise to a contribution to 
the potential that depends on the orientation of the dipoles and the inverse 
cube of the intermolecular separation. An appropriate potential for such 
molecules is the Stockmayer potential 



Here 8 is a function of the angular orientation of the molecules and also the 
strength of the dipole moment. It has not as yet been possible to calculate 
transport properties for this angle- dependent potential. However, collision 
integrals have been calculated for Stockmayer potentials modified by assuming 8 
a constant . These potentials are shown in Fig. 5; at the extremes of 8 = ±2.5 
the long-range part of the potential is completely dominated by the inverse cube 
term. (For comparison the Lennard-Jones (12-6) potential is shown as a dashed 
curve (8 = 0). ) 

These collision integrals may be applied to polar gases if it is assumed 
that the relative orientation of the colliding dipoleB remains fixed through 
the important part of the collision trajectory around the distance of closest 
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approach. The collision integrals are then averaged over all possible relative 

g 

orientations . Computed and experimental viscosities of steam (5 = 1.2) are 
shown in Fig. 6. The fit of the data is not perfect - the experimental data 
show a temperature dependence that is somewhat steeper that the theoretical 
prediction - but it is an improvement over the fit obtained by using a simple 
Lennard-Jones (12-6) potential. 

Collision integrals have also been computed for the range of Morse potentials 
shown in Fig. 7^. The potentials with the very broad wells can be applied to 
interactions between atoms corresponding to chemically bound molecules (for 
example, the interaction between two hydrogen atoms, which corresponds to 
the hydrogen molecule). The potentials with the narrow wells have been used to 
fit data on nonpolar, valence- saturated molecules. 

Collision integrals also have been computed for a number of other potentials 
including the repulsive exponential potential' 1 ' 1 (for molecules at high tempera- 
tures, as well as nonbonding interactions between atoms or free radicals), the 

12 

shielded coulombic potential (for repulsive interactions in ionized gases), 

several (12-6-4) potentials 15 (for ion-neutral interactions), several inverse 

power attractive and repulsive potentials 1 ^, and also estimates for the attractive 
15 , 

exponential (for low-temperature interactions between atoms or free radicals 
corresponding to bound molecular states). 

Thus the collision integrals are now available for a wide variety of 
potentials. Hence, if the details of a potential energy curve are known, a 
good match to it can be selected from the "library" of calculations already 


at hand. 
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HEAT CONDUCTION IN CHEMICALLY REACTING GASES 

At high temperatures piany gases are partially dissociated and undergo a 
variety of chemical reactions. In reacting gases, heat transport may he con- 
siderably larger than in "frozen" (nonreacting) mixtures. Large amounts of heat 
can be carried as chemical enthalpy of molecules that diffuse because of concen- 
tration gradients. These gradients exist, in turn, because the gas composition 
varies with temperature. For example, in a gas that absorbs heat by dissociating 
as the temperature is raised, heat is transported when a molecule dissociates in 
the high-temperature region and the fragments diffuse toward the cooler region. 

In the low-temperature region the fragments recombine and release the heat 
absorbed at high temperature. 

When chemical reaction rates are very high, chemical equilibrium can be 
assumed to exist locally throughout a gas mixture. It Is then possible, by 
differentiating the equilibrium relationships, to relate the concentration 
gradients to the temperature gradient. In this event one can define an equilib- 
rium thermal conductivity A g independent of apparatus geometry: 

A e = A-p + A r (10) 

where A^ is the conductivity in the absence of reaction (the "frozen" 
thermal conductivity) and A r is the augmentation due to the reactions. 

A general expression for the thermal conductivity due to chemical reactions 

1 f ? "17 

hap been developed ' that is applicable to mixtures involving any number of 
reactants, inert diluents, and chemical equilibria, provided chemical equilib- 
rium exists locally in the temperature gradient. For a simple dissociation of 
the type A ^ nB the thermal conductivity due to chemical reaction is 
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( 11 ) 
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Here D^-g is the "binary diffusion coefficient between components A and B, 

AH is the heat of reaction, and x^, Xg are the mole fractions of the com- 
ponents. Note that unless both species are present, A r is zero. Furthermore, 
since in a dissociating gas the gas composition varies with pressure, we expect 
the heat conductivity to vary with pressure also. This is in contrast to the 
behavior of nonreacting gases, for which the heat conductivity is independent 
of pressure. 

Experimental and theoretical conductivities for the NgO^ ■** 2N0g system 

18 

at one atmosphere are shown in Fig. 8 . The dashed curve indicates the frozen 

conductivity. Thus is the major contribution to the heat conductivity; at 

the maximum (where the mass fractions of Ng0 4 and NOg are equal) the con- 
ductivity is comparable to that of a light gas such as helium. 

The theoretical expression for a system involving two reactions has been 
19 

tested for the case of hydrogen fluoride vapor. At moderate pressures the 
PVT behavior of hydrogen fluoride can be described in terms of a monomer-hexamer 
equilibrium, while low pressure data suggest a dimer as well. Although the 
actual state of the vapor is uncertain, it appears that at low and moderate 
pressures the equilibria 


2HF 

* (hf) 2 

AHg = 7.4 kcal 

6HF 

- (HF)g 

AH g =40.5 kcal 


serve to specify the system rather well. 

20 

Computed and experimental thermal conductivities are compared in Fig. 9. 

The solid line was computed assuming both dimer and hexamer equilibria, whereas 
the dashed line was computed considering only the hexamer equilibrium. Note the 
extreme pressure dependence of the thermal conductivity. The maximum conductivity 
is more than three times that of hydrogen at the same temperature and some 33 
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times the frozen thermal conductivity expected in the absence of reaction. The 

inclusion of a dimer equilibrium markedly Improves the agreement between theory 

and experiment in the low-pressure region. 

The experimental studies on nitrogen tetroxide and hydrogen fluoride prove 

the validity of the theoretical expressions for thermal conductivity of reacting 

gases in chemical equilibrium. Recently the theory has been successfully applied 

21 

to data for the PCl^ PCI,. + Clg equilibrium . 

Thus far we have considered systems where the chemical reactions are so 

rapid that chemical equilibrium prevails locally at all points in the gas mixture. 

Let us now consider the reduction of heat transport caused by reduced reaction 

22 

rat.es. A general expression has been derived for the apparent "thermal con- 
ductivity" of reacting mixtures in which a single reaction proceeds at a finite 
rate. In contrast to systems where reaction rates are either very high or very 
low, it is found that heat conduction depends on the geometry and scale of the 
system and also the catalytic activity of the surfaces. 

For a plane parallel plate geometry, with one surface noncat alytic and the 
other surface a perfect catalyst, the effective "thermal conductivity" is 

A» « (12) 

A f + A r 522^-2 

where 
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Here is the chemical reaction rate at equilibrium (that is the total rate 
in either direction - not the net rate, which is zero, of course), and Z is 


the distance between the plates. For simple systems it can be shown that 
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If the diffusion time is short in comparison with the chemical relaxation 

time T c h em > 9 -*■ tanh cp/cp ■* 1, and X* -*• A^. In other words, the concentra- 
tion gradients are washed out by diffusion and the frozen conductivity is obtained. 
On the other hand if the chemical time is short, the concentration gradients are 
maintained, cp -* 1 , tanh cp/cp *♦ 0, and X* -* A 0 . 

18 

The theory has been applied to the low-pressure measurements on the 
Ng0 4 £ 2N0g system, as shown in Fig. 10, The upper and lower dashed curves are 
respectively, the computed equilibrium and frozen conductivities. The remaining 
two curves are calculations of X* for various reaction rates, assuming neg- 
ligible chemical reaction on the surfaces. Near atmospheric pressure the experi- 
mental data lie on the equilibrium conductivity curve, which is in agreement 
with Fig. 8. 

At low pressures the dissociation of Ng0 4 is a fast bimolecular reaction: 

N 2 0 4 + M - 2N0g + M 

The curve marked M = Ng has been calculated by using the rate data of Carrington 
and Davidson for the dissociation of N2O4 in nitrogen. The solid curve has 
been calculated by assuming the second order rate is sevenfold greater when 
undiluted N 2 0 4 - NOg mixtures dissociate* this is in conformity with the 
experiments of Bauer and Gustavson . The agreement between theory and experi- 
ment is very satisfactory. 

THERMAL CONDUCTIVITY OF POLYATOMIC GASES 

It is convenient to discuss the thermal conductivity of a polyatomic gas 
in terms of its relationship to the viscosity through the dimensionless ratio 

f = AM/qC^. (14) 

Here M is the molecular weight and Cy is the constant volume molar heat 
capacity. According to ultras Amplified kinetic theory, f = lj however, the 
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rigorous Chapman- Enskog theory for monatomic gases predicts that f should he 
very nearly 5/2. This is due to the fact that translational energy is a function 
of molecular velocity; the molecules possessing the most energy are the most 
rapid, have the longest mean free paths, and hence make an enhanced contribution 
to the heat transport. Indeed, experiment confirms that f is about 2.5 for 
the noble gases. This is illustrated in Fig. 11 where data for argon are shown 
over a temperature from about 100° to 300° K. The data points represent measure- , 
ments^ 5 that provides a direct determination of 

For polyatomic gases, f is less than 2.5 and tends to be smallest when the 
molar heat capacity is largest and originates mostly from the internal energy 
modes. Consequently, Eueken^ 7 suggested that the transport of translation and 
internal energy be considered separately, and proposed 

fC v = f-trans ^v, brans + ^int ^int (15) 

(C J and C. . are the translational and internal contributions to the 

v v,trans int 

total heat capacity C T . ) Eueken assumed f-fcrans = ®/2, by analogy with the 
monatomic gases. However, because there is little correlation between molecular 
velocity and internal energy, Eueken assumed f = 1 (the result of the ultra- 
simple theory that neglects the velocity-translational energy correlation). 

OO 

Ubbelohde pointed out that molecules with excited internal energy states 
may be regarded as different chemical species and that the flow of internal 
energy can be considered as energy transport due to diffusion of the excited 
states. This concept leads to the result f = pD/q, so that 

# " ' r ' " ' ir ,r ' _ ' ' 

The method involves measurement of the adiabatic recovery temperature T r 

attained on a flat plate in a high-velocity subsonic gas stream. The recovery 
temperature is related to the stream temperature T„ and total temperature 
T-k through the recovery factor r = (T r -T g )/(T^-T s ). To a very good approx- 
imation f = Y"/ r ^ where y is the specific heat ratio. 


For many realistic force laws pD/rj « 1.3 over a large temperature range. 

To justify this modified Eucken approximation [Eq. (16)] it is tacitly assumed 
that inelastic collisions are rare. This is necessary in order that the 
translational velocity distribution function should not be unduly perturbed, 
so that the translational conductivity can be related to the viscosity as in 
the case of the noble gases. On the other hand, there must be enough inelastic 
collisions to maintain the internal energy states in equilibrium with the 
local temperature. 

29 

Mason and Monchick have recently derived explicit expressions for 
ft rans an< ^ fint ^ rom "the formal kinetic theory of polyatomic gases. By 
systematically including terms involving inelastic collisions they obtained 
the modified Eucken expression as a first approximation, and, as a second 
approximation, an expression dependent on the relaxation times for the various 
internal degrees of freedom. For nonpolar gases their result may be written 


2 C 


f MM C v = ¥ R + C int ' I ( I " rf ) 


(17) 


Here C rafc is the rotational contribution to the heat capacity (R for linear 
molecules, ^ R for nonlinear molecules) and Z rQ ^. is a collision number for 
rotational relaxation: 


^rot = T rot / l coll 


-,/ T r 


4 RT rot 


(18) 


where t ^ is the rotational relaxation time and T co qi = ( Jt / 4 )( T 1 AO is the 
mean time between collisions. 

25 . ,. 30, 31 _ , , 25 

Experimental f values for nitrogen , carbon dioxide , and hydrogen 


are shown in Figs. 12, 13, and 14. The upper dashed curves correspond to the 
modified Eucken approximation [Eq. (16)], while the solid lines have been 
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calculated from Mason and Monchick’s expression [Eq* (17)] by assuming the 
temperature- independent collision numbers shown in the figures. The Eucken 
approximations (f-trans = 2,5, ?±nt = are as well for completeness. 

The data on nitrogen and carbon dioxide lie midway between the Eucken and 
modified Eucken approximations. In contrast, the hydrogen data scatter about 
the modified Eucken approximation (Z = co). Indeed, the hydrogen molecule is 
unique in that exchange between translational and rotational energy is not 
easy; collision numbers of a few hundred are predicted from theory and observed 
experimentally. 

Collision numbers determined from f, or recovery factor, are compared with 
other measurements in Table I. The values from f are generally in accord 
with those obtained by the other techniques, within the admittedly rather large 
uncertainties associated with such determinations. The only serious disagree- 
ment is in the case of carbon dioxide, where acoustical measurements indicate 
that 16 collisions are required for relaxation. It seems likely that the 
acoustic result is in error; it is difficult to see any theoretical reason why 
carbon dioxide should relax so slowly. 

A classical theory for the rotational relaxation of molecules with attrac- 
tive intermolecular forces (rough spheres and spherocylinders surrounded by 

70 

square wells) has recently been developed by Sather and Dahler . In the case 
of rough spheres the rotational relaxation time is 


-1 

rot 
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(4l/mo 2 ) 


/ itkT \ 

\ m ) 


g(cr) 


[l + (4l/mo^) ]' 

Here n is the number of molecules per cm^, I is the moment of inertia, m 
is the mass, a is the diameter of the rough sphere core, while g(o) is the 


(19) 


value of the radial distribution function at a. In the low-density limit, 
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g ( cf ) = exp (e/kT), where e is the depth of the well at a. 

Equation (19) possesses two noteworthy features. First, since the 
quantity 4l/mcr^ is generally less than 0.2, the relaxation time depends only 
weakly on a, the position of the repulsive core. Secondly, since there is no 
rotational energy transfer accompanying the velocity impulse at the outer edge 
of the potential well, the only contribution to the rotational relaxation stems 
from the impulse at the core. Consequently, the width of the well is unimpor- 
tant; in fact, Eq. (19) should apply to a Sutherland potential with a rough 
core - that is, a potential with an inverse sixth power attractive portion to 

f 

account for Van der Waals forces. We might hope that this model would be suit- 
able for molecules that are approximately spherical. 

After combining Eqs. (l), (18), and (19) we find 


Z -1 ~ 5£i 
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exp( e/kT) 


( 20 ) 
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(The term 4l/ma^ in the denominator of Eq. (19) has been neglected. ) Some 
experimental results are compared with predictions of Eq. (20) in Table II. 
Methane, carbon tetrafluoride, and sulfur hexafluoride are approximately 
spherical, as shown in Fig. 15; for these molecules the agreement between theory 
and experiment seems very good indeed. The mass distribution parameter 
4I/mcrfr ’ ' varies more than threefold between methane and carbon tetra- 
fluoride; thus the large collision number of methane is probably a direct 
consequence of the molecule's small moment of inertia. Note that the effect 
of attractive forces [exp (e/kT)] is appreciable and roughly doubles the 
transition probabilities. 

The calculations for ethylene and ethane were carried out using average 

moments of inertia. The Z~\ values calculated for ethane are in close accord 

rot 
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with experiment, hut in the case of ethylene the agreement is not so good. 

From Fig. 15 it is apparent that the ethane molecule is approximately spherical, 

whereas the ethylene structure is definitely less compact and symmetric. 

It appears, then, that the classical kinetic theory for a rough spherical 

molecule with attractive forces may provide a lower limit to the collision 

probability for rotational relaxation, Z~\, for nonlinear molecules. Molecules 

rot 

such as CH4, CF4, SFg, and 03% are reasonably represented by this model, whereas 

less symmetric molecules such as CgH^ have large transition probabilities, or 

31 33 

shorter relaxation times. As a matter of fact the theory of Brout ’ for the 
relaxation of diatomic molecules such as nitrogen and oxygen also indicates 
that the deviation of the intermolecular potential from spherical symmetry is 
an important parameter. 

Thus we conclude that Mason and Monchick’s approximate theory for the heat 
conductivity of polyatomic gases is at least qualitatively correct. Furthermore, 
the following factors seem of profound importance in determining rotational 
relaxation times for nonpolar gases: 

(1) The mass distribution (characterized by 4l/ma^) 

(2) The strength of the intermolecular attractive forces (characterized 

by e/kT) 

(3) The deviation of the molecular force field from spherical symmetry 

CONCLUDING REMARKS 

In the past decade there have been substantial advances in methods for 
calculating the transport properties of gases and gas mixtures. Collision 
integrals have been calculated for a considerable variety of realistic potential 
energy functions. The effects of chemical reaction and chemical rate phenomena 



16 


seem to "be well understood, and there is gratifying accord between theory and 

29 

experiment. Finally, the theory of Mason and Monchick shows great promise 
as a description of the behavior of polyatomic gases. 

What remains for the future? First, it is to be hoped that the relaxation 
theory of Mason and Monchick will be extended to gas mixtures, and that there 
will be further testing of the theory, both for pure gases and gas mixtures. 
There is room for more work on the heat conductivity of polar gases. The con- 
ductivities of highly polar gases seem anomalously low in relation to their 

OQ 

viscosity (in other words, f values are small), and it has been suggested 

that this effect is largely due to a resonant exchange of rotational quanta, 

presumed probable on grazing self -collisions of polar molecules. However, 

34 

results of experiments designed to test this notion have been somewhat 
ambiguous. Finally, theoretical studies on the kinetic theory of nonspherical 
molecules should be encouraged. 
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TABLE I. - COMPARISON OF COLLISION NUMBERS FOR ROTATIONAL RELAXATION 


Gas 

f 

(determined 
from 
recovery- 
factor ) 

Acoust ical 

Low- 

pressure 

thermal 

conduc- 

tivity 

Shock 

thickness 

Impact 

tube 

% 

7 

5, 4-6 


5.5 

<7, <14 

°2 

12 

2-4, 4.1, 12, 
14,12-30 

20 

7 


H 2 

Large 

240-360 

300 

>150 

160, 310 

co 2 

2.4 

16 




ch 4 

9 

14-17 





TABLE II. - EXPERIMENTAL DATA FOR SOME NONLINEAR MOLECULES 


COMPARED WITH ROUGH SPHERE THEORY 



Methane 

Carbon 

tetrafluoride 

Sulfur 

hexafluoride 

Ethane 

Ethylene 

4l/,»0 2 n< 2 ’ 2 >* 

0.05 

0.16 

0.14 

0.09 

0.08 

exp( e/kT) 

1.7 

1.6 

2.1 

2.1 

2.2 

Z - -'- calc 

.11 

.33 

. 38 

.25 

.23 

Z" 1 expt 

.11 

.33 

.36 

.25 

.42 

Z expt 

9 

3.0 

2.8 

4 

2.4 
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Fig. 5. - Modified Stockmayer potentials for which collision integrals 
have been computed. 
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Fig. 7. - Morse potentials for which collision integrals have been computed. 
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Fig. 8. - Thermal conductivity of the N 2 0 4 -N0 2 system at atmospheric 
pressure. 
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Fig. 9. - Thermal conductivity of hydrogen fluoride vapor 
at 267. 7° K. 
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Fig. 15. - Configurations of nonlinear molecules. 
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